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Abstract — Power control is a fundamental task accomplished 
in any wireless cellular network; its aim is to set the transmit 
power of any mobile terminal, so that each user is able to 
achieve its own target SINR. While conventional power control 
algorithms require knowledge of a number of parameters of 
the signal of interest and of the multiaccess interference, in 
this paper it is shown that in a large CDMA system much of 
this information can be dispensed with, and effective distributed 
power control algorithms may be implemented with very little 
information on the user of interest. An uplink CDMA system 
subject to flat fading is considered with a focus on the cases in 
which a linear MMSE receiver and a non-linear MMSE serial 
interference cancellation receiver are adopted; for the latter 
case new formulas are also given for the system SINR in the 
large system asymptote. Experimental results show an excellent 
agreement between the performance and the power profile of the 
proposed distributed algorithms and that of conventional ones 
that require much greater prior knowledge. 

I. Introduction 

In multiuser wireless communication systems, mobile users 
vary their transmit powers so as to counteract channel gain 
variations and achieve their target Signal-to-Interference plus 
Noise Ratios (SINRs); this task is usually referred to as power 
control [1]. In code division multiple access (CDMA) wireless 
systems, power control, possibly coupled with the use of 
multiuser detectors, is thus used in the uplink to combat the 
near-far effect, to manage interference levels, and to minimize 
the overall power radiated by the system. A considerable 
amount of work has been done on power control algorithms 
for cellular networks, see, e.g. [l]-[3], to cite a few; these 
algorithms, based on measurements taken at both the mobile 
station (MS) and at the base station (BS), give as output the 
transmit power for each terminal. In power control procedures, 
estimates of several parameters such as the channel gains and 
the SINR for each user are usually needed, or, alternatively, 
recursive algorithms are adopted (see, e.g. [3]), which however 
are affected by slow convergence speed and excess steady-state 
error 

In recent years a new mathematical tool has emerged in 
the analysis of CDMA systems, i.e. the so-called large system 
analysis, first introduced in [4]. In short, [4] has revealed that, 
in a CDMA system whose processing gain and number of 
users both increase without bound with their ratio fixed, and 
with randomly chosen unit-norm spreading codes, the SINR 
at the output of a linear minimum mean square error (MMSE) 



receiver converges in probability to a non-random constant. In 
particular, denoting by K the number of active users, by TV the 
system processing gain, by the additive thermal noise power 
spectral density (PSD) level, and by Ep[-] the expectation with 
respect to the limiting empirical distribution F of the received 
powers of the interferers, the SINR of the MMSE receiver for 
the k-th user, say jk, converges in probability as K,N 
oo, K/N = a = constant, to 7^ the unique solution of the 
equation 
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with Pk the received power for the fc-th user. Interestingly, 
the limiting SINR depends only on the limiting empirical 
distribution of the received powers of the interferers, the load 
a, the thermal noise level and the transmit power of the user 
of interest, while being independent of the actual realization 
of the received powers of the interferers and of the spreading 
codes of the active users. Large system analysis is now a well- 
established mathematical tool for the design and analysis of 
communication systems (see, e.g., [5], [6], to cite a few). 

In this paper, we show how large system analysis can 
be used to design distributed power control algorithms that 
need very little prior information (i.e. the channel gain for 
the user of interest) to be implemented, for both the cases 
in which linear MMSE detection and non-linear interference 
cancellation MMSE detection are used at the receiver The 
contributions of this paper can thus be summarized as follows. 

- Assuming that linear MMSE detection is used at the 
receiver, we propose a new distributed power control 
algorithm requiring little prior information and based on 
large system analysis. 

- Extending the approach in [4], we give an expression for 
the limiting SINR in the case in which non-linear serial 
interference cancellation MMSE (SIC/MMSE) detection 
is adopted at the receiver 

- Assuming that SIC/MMSE detection is used at the re- 
ceiver, we propose a new distributed power control algo- 
rithm based on large system approximations. 

- We show how the proposed algorithms can be included 
in utility-maximizing non-cooperative games in order to 
achieve energy-efficiency in wireless data networks [7]. 



Our numerical results will show that the proposed algorithms 
are effective and are able to approach very closely the power 
profile predicted by much more complex and information- 
demanding power control algorithms. 

II. System model and problem statement 

Consider the uplink of a K-usei synchronous, single-cell, 
direct-sequence code division multiple access (DS/CDMA) 
system with processing gain N and subject to flat fading. After 
chip-matched filtering and sampling at the chip-rate, the A^- 
dimensional received data vector, say r, corresponding to one 
symbol interval, can be written as 
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wherein pk is the transmit power of the fc-th useiQ, bk E 
{—1,1} is the information symbol of the fc-th user, and hk 
is the reajl channel gain between the fc-th user transmitter and 
the base station; the actual value of hk depends on both the 
distance of the fc-th user's mobile from the base station and 
the channel fading fluctuations. The A^-dimensional vector Sk 
is the spreading code of the fc-th user; we assume that the 
entries of Sk are binary-valued and that sJ^Sk = WskW^ — 1, 
with (•)^ denoting transpose. Finally, n is the thermal noise 
vector, which we assume to be a zero-mean white Gaussian 
random process with covariance matrix cf'^In, with Im the 
identity matrix of order N . 

Given the system model (|2]l, a number of strategies are avail- 
able to detect the data symbols hi, ... ^bK', in the following 
we briefly review the hnear MMSE receiver and the non-linear 
SIC/MMSE receiver. 

A. Linear MMSE detection 

Consider the case in which a linear receiver is used to 
detect the data symbol bk, according, i.e., to the decision 



rule bk 



sign 



with bk the estimate of bk and 
dk the A^-dimensional vector representing the receive fil- 
ter for the user fc. Then, it is easily seen that the lin- 
ear MMSE receiver is the one corresponding to the choice 

dk = ^/Pkhk i^SHPH'^S'^ + (t'^In^ Sk , wherein S = 
[si, . . . , Sk] is the A^ x iC-dimensional spreading-code matrix, 
and P and H are K x A'-dimensional diagonal matrices, 
whose diagonals are [pi,...,pk] and [hi,...,hK], respec- 
tively. For linear detectors it is also meaningful to define the 
output SINR, which, for the fc-th user is written as 



Ik 



Pkhlidlsk 



^p,hl{dls,f 



(3) 



'To simplify subsequent notation, we assume tliat the transmitted power 
Pk subsumes also the gain of the transmit and receive antennas. 

-For simplicity, we assume a real channel model; however, generalization 
to practical channels with I and Q components is straightforward. 



B. Non-linear SIC/MMSE detection 

Consider now the case in which non-linear decision feed- 
back detection is used at the receiver We assume that the 
users are indexed according to a non-increasing sorting of 
their channel gains, i.e. we assume that hi > h2 > . . . , h^- 
We consider a serial interference cancellation (SIC) receiver 
wherein detection of the symbol from the fc-th user is made 



dlvk 



wherein 



according to the rule bk — sign 

Otherwise stated, when detecting a certain symbol, the 
contributions from the data symbols that have already been 
detected are subtracted from the received data. The output 
SINR for user fc, under the assumption of correctness of past 
decisions, is now written as 
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Upon defining Sk = [sk, ■ • • , sk], Pk = diag(pfe, . . . ,pk) 
and Hk — diag{hk, . . . , hx), it is easy to show that 
SIC/MMSE detection corresponds to the choice dk ~ 
JpihkiSkHkPkHlSl + a^lN^^Sk . 



C. Problem statement 

Given the data model dU, we are interested in the following 
problem: find the transmit power pk G [0, Pmax] for each 
user fc, so that the SINR jk equals a given target value 7, 
with Pmax the maximum power that each user in the system 
is allowed to transmit. Note that this problem finds numer- 
ous applications. As an example, in circuit-switched wireless 
cellular voice communications, wherein the primary goal is 
signal intelligibility, the SINR is required to be always above 
a given intelligibility threshold [1]. In wireless packet-switched 
networks, instead, the primary goal may be to maximize the 
system throughput, or, if battery-life of the mobile terminals 
is a dominant issue, the system throughput for each unit or 
energy drained from the battery [2], [7]. In all cases, however, 
it can be shown that this goal translates into the requirement 
that each user's SINR equals at least a certain value. 

In the sequel, we show how large system analysis leads 
to power control algorithms that may be implemented in a 
distributed fashion and that require knowledge of the channel 
for the user of interest only. 

III. Power Control for linear MMSE detection 

As an introductory step in our algorithm, we begin by 
illustrating a simple power control algorithm derived from 
[4]. We have seen that in a large CDMA system the fc- 
th user's SINR converges in probability to the solution to 
Eq. ([T]l. Heuristically, this means that in a large system, and 
embracing the notation of the previous section, the SINR 7fc 
is deterministic and approximately satisfies 

hlPk 
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Now, as noted in [4], if all the users must achieve the same 
common target SINR 7, it is reasonable to assume that they 
are to be received with the same power, i.e. the condition 
hfpi = ft,2P2 — ■ ■ ■ h\pk — Pr, , is to be fulfilled. 
Substituting the above constraint in (|5]l and equating Q to 
7 it is straightforward to come up with the following relation 
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7(7 



1 



1 



Pk 



"fa 



1+7 



(6) 



1+7 



wherein, we recall, a = K/N, and the relation a < 1 + 
1/7 must hold. Eq. (|6]l, which descends from eq. (16) in [4], 
gives a simple power control algorithm that permits setting 
the transmitted power for each user based on the knowledge 
of the channel gain for the user of interest only. The above 
algorithm, however, does not take into account the situation 
in which, due to fading and path losses, some users end up 
transmitting at their maximum power without achieving the 
target SINR, and indeed our numerical results to be shown in 
the sequel will prove the inability of (|6]l to predict with good 
accuracy the actual power profile for the active users. 

In order to circumvent this drawback, we first recall that in 
[8] (see also [9]) the following result has been shown: 
Lemma: Denote by F{-) the cumulative distribution func- 
tion ( CDF) of the squared fading coefficients hf, and by 



, , h"^^-^ the vector of the users' squared fading 



coefficients sorted in non-increasing order. Then we have that 



K^p^ converges in probability (as K 
y£^l,...,K. 



00) to F-i 



The above lemma states that if we sort a large number 
of identically distributed random variates, we obtain a vector 
that is approximately equal to the uniformly sampled version 
of the inverse of the common CDF of the random variates. 
Accordingly, in a large CDMA system each user may indi- 
vidually build a rough estimate of the fading coefficients in 
the network and thus will be able to predict the number of 
users, say U2, that possibly will wind up transmitting at the 
maximum power. Indeed, since, according to ^ each user is 
to be received with a power Pfj, the estimate U2 of the number 
of users transmitting at the maximum power is given by 
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with m( ) the step-function. It is also natural to assume that the 
users transmitting at P,nax will be the ones with the smallest 
channel coefficients, i.e. the squared channel gains of the users 
transmitting at the maximum power are well approximated 
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1, 



by the samples F^^ {^^)^ with £ = K 
As a consequence, the generic k-th user will be affected by 
Ml = K — U2 users that are received with power Pr (these are 
the Ml users with the strongest channel gains and that are able 
to achieve the target SINR 7), and by U2 users that are received 
with power PmaxF^^ ("^)' with £ ^ K - U2 + I, ■ ■ ■ , K. 
Denoting by Pk the received power for the fc-th user, Eq. ^ 



can be now written as 
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Now, assuming for the moment that user k is able to achieve 
its target SINR, i.e. that Pk — Pr, we can make the 
approximation p^^p^^^ ~ ' whereby equating (O to 

the target SINR 7 we have 



Pk 



ui Pk 

N 1+7 



1 X^'^ 



K PkP.-...F-^{^) 
t=K-U2 + l Pfc+P„„,F-i(i^^)7 



7 ■ 



(9) 

The above relation can now be solved numerically in order to 
determine the receive power Pk for the fc-th useo Finally, the 
actual transmit power for the fc-th user is set according to the 
rule 

Pk=min{Pk/hl, P„,,^} . (10) 

The proposed algorithm may be summarized as follows. 
First, the number of users transmitting at the maximum power 
is estimated according to (|7]i. Then, the desired receive power 
for each user is computed by solving (|9]l. Finally, the transmit 
power for the fc-th user is determined according to relation 
( [Tol l. Note that this algorithm requires knowledge only of the 
channel gain for the user of interest. 

IV. Power control for SIC/MMSE detection 

Let us now consider the case in which non-linear 
SIC/MMSE detection is used at the receiver, and let us thus 
assume that users are ordered according to a non-increasing 
power profile. In this case the following theorem can be 
proved. 

Theorem: Let 7^ be the (random) SINR of the SIC/MMSE 
receiver for the k-th user; let Pk be the received power for the 
k-th user, and assume that previously detected symbols have 
been perfectly cancelled. As K,N ^ +00, with K/N = a, 
7fc converges in probability to 7^, the unique solution of the 
equation 
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with Oik — [K ^ k + l)/N and Ep\p^p^ [■] denoting expecta- 
tion with respect to the empirical distribution of the received 
powers not larger than Pk. 

Proof: The proof is omitted for the sake of brevity. 

To corroborate the statement of the above theorem, in Fig. 
1 we report the asymptotic SINR and 100 actual realizations 
of the SINR corresponding to random realizations of the 
spreading codes and of the received powers, assumed to follow 
a Rayleigh distribution. A system with processing gain ~ 
256 has been considered here, and it is seen that the random 
SINR reaUzations are spread around their asymptotic value. 

Actually this equation gives the desired receive power for each user, and 
thus in a centralized power control algorithm needs to be solved only once. 



Using the notation of Section II, a heuristic reformulation 
of the above theorem states that in a large system the fc-th 
user's SINR 7^ is deterministic and approximately satisfies 
the equation 
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Now, the above expression can be used to derive a simple 
and effective power control algorithm. Let us consider the K- 
th user first, i.e. the one with the smallest channel gain. The 
information symbol from this user will be the last one to be 
detected, thus implying that, under the assumption of error-free 
detection of previous bits, and denoting by Px the received 
power for this user, the relation P^ /cr^ = 7 should hold. As 
a consequence, the transmit power for the K-\h user is set to 

Pk = min{7(T^//l|-,P,nax}- 

Consider now user K—1; on denoting by Pk-i its received 
power, we have that 

Pk-1 
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" ' N PK-l+PKh^ 

Now, we may reasonably assume that Hk and Hk-i are 
approximately equal (recall that these are the two smallest 
channel gains), and thus that Pk-i ~ Pk- Moreover, in a 
distributed approach, we can substitute Kk with its estimate 
fiK given by the Lemma of the previous section. As a 
consequence, ( fT3] l can be approximated as 

Pk-1 
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^ Pk+PkIi^^ 

After solving the above equation for Pk-i, we can set the 
transmit power of the {K ~ l)-th user according to pk-i ~ 

min{Px_l//l|._;^, Pmax}- 

In general, for the generic fc-th user, we have that the 
following relation must hold 

Pk 
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Since it is reasonable to assume that Pk ~ Pk+i, and replacing 
the channel gains with their estimates, the above equation can 
be re-written as 
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Solving with respect to Pk we have 
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and the transmit power for the fc-th user is 

Pk = mm{Pk/hl , P,„ax} ■ 



(16) 



(17) 



(18) 



In summary, the proposed algorithm proceeds as follows. 
For centralized implementation, equations dTTI l and dTsl l are 




Fig. 1. Randomly generated SIC/MMSE SINR compared to the asymptotic 
limit, for the case of exponentially distributed powers (Rayleigh fading). 




Fig. 2. Transmitted power profile across users for the proposed algorithm, 
the conventional power control algorithm [1] and the profile derived by the 
algorithm in [4], for both cases of linear MMSE detection and SIC/MMSE 
detection. 



sequentially implemented for k = K, K —1, . . . ,2,1. Alterna- 
tively, for decentralized implementation, the generic j-th user 
must compute equations ( fTTj l and the equation 



Pk = 
K,K - 1, 



inm{Pk/hl , Pinax} , 



(19) 



for fc — K, K — 1, . . . , j + 1. Finally, equations ([TtT i and 
(fTSl l are implemented for k = j. Note that in the distributed 
implementation each user needs to know its own channel gain 
and also its order in the data detection sequence at the receiver, 
i.e. it must know how many users interfere with it at the 
receiver. However, no information on its uplink SINR or on 
the parameters of the multiaccess interference is needed. 

V. Numerical Results 

We consider an uplink DS/CDMA system with processing 
gain N = 128, and assume that the packet length is M — 120, 
and that the target SINR is 7 = 6.689 =8.25 dB. According 
to [7], it can be seen that achieving the target SINR 7 in the 
considered scenario leads to a maximization of the utility, i.e. 
of the ratio between the packet success rate and the transmit 
power 

A single-cell system is considered, wherein users may have 
random positions with a distance from the BS ranging from 
10m to 1000m. The channel coefficient hk for the generic fc-th 
user is assumed to be Rayleigh distributed with mean equal to 
d^^, with dk being the distance of user fc from the BS. As to 
the thermal noise level, we take — 2 - lO^^W/Hz, while the 
maximum allowed power Pk.max is — 25dBW. Fig. 2 reports 



the transmitted power profile across users for the algorithm 
proposed in Section III and IV, for the algorithm derived by 
Eq. (16) in [4] (i.e. eq. (|6]l), and for a conventional algorithm 
(see [1]) that is non-adaptive and requires a substantial amount 
of prior information. It is seen that the proposed algorithms 
are capable of reproducing the optimal power profile with 
very good accuracy, while, on the contrary, the algorithm 
descending from paper [4] overestimates the required transmit 
powers and does not exhibit good performance. It is seen 
clearly that adopting a SIC/MMSE receiver yields considerable 
savings in transmit power needed to achieve a certain target 
SINR. 

While Fig. 2 shows the result of only one simulation trial 
(note however that similar behavior has been observed in every 
case we considered), the remaining three figures report results 
coming from an average over 1000 independent trials. Figs. 3 
- 5 show the achieved average utility (measured in bits/Joule), 
the average per-user transmit power and the average achieved 
SINR at the receiver output for the conventional power control 
algorithms (for both linear MMSE and SIC/MMSE detection) 
[1], for the proposed algorithms (for both linear MMSE and 
SIC/MMSE detection), and for the power control algorithm 
derived by paper [4]. These results show that the proposed 
algorithms achieve performance levels practically indistin- 
guishable from those of the standard algorithms, while the 
algorithm ^ achieves a much smaller utility. From Fig. 5 it 
is however seen that the algorithm ^ achieves an output SINR 
larger than that of any other algorithm considered: this should 
not be interpreted as a sign of good performance. Indeed, in 
the considered scenario the aim of the power control algorithm 
is to make each user operate at a SINR equal to 7. 

VI. Conclusions 

This paper has considered the design of distributed power 
control algorithms for cellular CDMA systems based on 
asymptotic analysis, for situations in which either linear 
MMSE detection or a non-linear SIC/MMSE detection are 
used by the receiver. For the latter case, closed-form formulas 
for the limiting system SINR have also been developed. 

Overall, the proposed solutions achieve satisfactory perfor- 
mance, and the proposed approach is quite promising. Among 
the authors' current research efforts in this area is the extension 
of the proposed algorithms to the situation in which the 
received signals have been affected by multipath. 
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